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I.  INTRODUCTION 


A  solution  formulation  and  some  numerical  results  are  presented  for  beam 
motions  subjected  to  moving  loads.  Most  of  the  work  on  this  problem  has  been 
related  to  rail  and  bridge  design  (see,  for  example,  reference  1  and  many 
papers  cited  there  from  1910  to  1971).  However,  the  application  of  the  analy¬ 
sis  can  obviously  be  extended  to  tracks  for  rocket  firing  and  to  gun 
dynamics. ^ 

In  Section  II  of  this  report,  a  variational  formulation  for  a  moving 
force  problem  is  described.  Also  given  are  the  procedures  which  lead  to 
finite  element  matrix  equation.  A  detailed  description  of  the  treatment  of  a 
concentrated  moving  force  is  given  in  Section  III.  The  variational  problem 
associated  with  a  gun  tube  dynamics  is  presented  in  Section  IV.  This  gun-tube 
problem  contains  the  moving  mass  problem  as  a  special  case.  Finite  element 
solution  can  be  derived  from  this  formulation,  but  the  details  of  this  more 
complicated  problem  is  omitted  from  the  present  report.  Some  of  the  numerical 
problem  are  reported  in  the  last  section  and  are  compared  with  results 
obtained  from  series  solutions. 


^Fryba,  L. ,  Vibrations  of  Solids  and  Structures  Under  Moving  Loads,  Noordhoff 
International  Publishing  Company,  Groningen,  1971. 

^Wu,  J.  J.,  "The  Initial  Boundary  Value  of  Gun  Dynamics  Solved  by  Finite 
Element  Unconstrained  Variational  Formulations,”  Innovative  Numerical 
Analysis  For  the  Applied  Engineering  Science,  R.  P.  Shaw,  et  al.,  Editors, 
University  Press  of  Virginia,  Charlottesville,  1980,  pp.  733-741. 
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II.  SOLUTION  FORMULATION  FOR  A  MOVING  FORCE  PROBLEM 


In  this  section,  the  solution  formulation  will  be  described  in  detail  for 
a  moving  force  problem.  The  moving  mass  problem  will  be  included  as  a  special 
case  of  a  more  general  problem  of  gun  motions  analysis  given  in  a  later 
section . 

Consider  a  vertical  force  P  moving  on  an  Euler-Bernoulli  beam.  The 
differential  equation  is  given  by 

Ely” ,f  +  pAy  =  PS(x-x)  (1) 

where  y(x,t)  denotes  the  beam  deflection  as  a  function  of  spatial  coordinate  x 
and  time  t.  E,  I,  A,  p  denote  elastic  modulus,  second  moment  of  inertia,  area 
and  material  density  respectively.  A  Dirac  function  is  denoted  by  <S,  x  =  x(t) 
is  the  location  of  P,  a  prime  (f)  denotes  differentiation  with  respect  to  x 
and  a  dot  (•),  differentiation  with  respect  to  t. 

Introducing  nondiraensional  quant  it ies 

AAA 

y  -  y/i  ,  x  -  x/i  ,  t  =  t/T  ,  (2) 

where  l  is  the  length  of  the  beam  and  T  is  a  finite  time,  within  0  <  t  < 

T,  the  problem  is  of  interest,  Eq .  (1)  can  be  written  as 

y""  +  y2y  =  Q<S(x-x)  (3) 


2 


The  hats  (A )  have  been  omitted  in  Eq.  (3)  and 


with 


(4) 


Boundary  conditions  associated  with  Eqs.  (1)  or  (2)  will  now  be  introduced  in 
conjunction  of  a  variational  problem.  Consider 

61  =  0  (5a) 


with 

11  *  *  —  “ 

I  =  f  j  [y"y*”  -  y2yy*  -  Q6(x-x)]dxdt 
0  0 

+  J  dt{k1y(0,t)y*(0,t)  +  l^y'  (0,t)y*'  (0,t) } 

+  Y2/1  dx{k5[y(x,0)  -  Y(x) ] y*(x, 1) }  (5b) 

0 


where  y*(x,t)  is  the  adjoint  variable  of  y(x,t).  If  one  takes  the  first 
variation  of  I  considering  y(x,t)  to  be  fixed: 

(fil)fiy-O  3  0  (5a,) 

and  consider  6y*  to  be  completely  arbitrary,  it  is  easy  to  see  that  Eq.  (5) 
is  equivalent  to  the  differential  equation  (3)  and  the  following  boundary  and 
initial  conditions. 
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y'M  (0 ,t)  +  kiy(0,t)  =  0 
y" (0,t)  -  k.2y'  (0,t)  =*  0 
y”' (l,t)  -  k3y(l,t)  =  0 

y"  (1  ,t)  +  k^y ' (1 ,t)  =  0 


a  <  t  <  1 


(6a) 


y(x.O)  =  0 


and 


0  <  x  <  1 


(6b) 


y(x,l)  -  k5[y(x,0)  -  Y(x)]  =  0 


Taking  appropriate  values  for  k^,  k2,  k3,  and  k4,  problems  with  a  wide  range 
of  boundary  conditions  can  be  realized.  The  initial  conditions  in  Eqs.  (6b) 
are  that  the  beam  has  zero  initial  velocity,  and,  if  one  takes  k5  to  be  °°  (or 
larger  number  compared  with  unity), 

y(x ,0)  =  Y(x) 

The  meaning  for  cases  where  k5  is  not  so,  need  not  be  our  concern  here. 

To  derive  the  finite  element  matrix  equations,  one  begins  with  Eq. 

(5a')  and  write 


(fiDfiy-0  =  0 


(7a) 


+  k3y(l  ,t)S.y*(l  ,t)  +  k4y'  ( 1  ,t ) Sy* '  ( 1  ,t )  ] 


(7b) 
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Introducing  element  local  variables 


or 


(i) 

5=5  =  Kx-i+1 

(8a) 

(i) 

n  =  n  =  Lt-j+l 


1 

X  =  -  (5+i-l) 

K 

(8b) 

1 

t  =  -  (5+j-l) 

L 


where  K  is  the  number  of  divisions  in  x  and  L,  in  t .  (A  typical  grid  scheme 
is  shown  in  Figure  2).  Equation  (7b)  can  now  be  written  as 

l  l  Jl  Il[~  y"(ij)6y*"(ij)  -  ---  y(ij)6y*(ij)]d5dn 

i=l  j  =  l  U  U  L  K. 

L  I  ki  K2 

+  I  /  dn  [—  y(ij)(0,n)6y*(1j)(0,n)  +  k2  --  y' (ij)(0,n)6y*' (1j)(0,n) 
j_l  0  L  L 

+  l  J1  “  [Y2k5(y(ij) (5 ,0)6y*(i j) (5 ,1))] 

i=l  0  K 


K  L  Q  !  !  _  _ 

*  I  I  -  /  /  6(x-x)6y*(1j)(5,n)d5dn 

i=l  j=l  L  0  0 

K  Y2k5  i 

+  l  - I  d5  [Y(1) (5)6y*(iL)(5 ,1)1  (9) 

i  =  l  K  0 


The  shape  function  vector  is  now  introduced.  Let 

y(ij)(5,n)  -  aT(5,n)Y(1j) 

y*(ij)(5,n)  =  aT(5 >h)Y*(ij)  =  Y*T(ij)a(5 ,h)  (10) 
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Equation  (9)  then  becomes 


,KJ 


Y2L 


I  I  «?*T(lj)  P-  *  '  B|  Y(1J) 
1*1  1*1  L  K  ~  ~ 


L  kl  k2K2 

+  l  «y*T(ij)  (--  Bl  +  —  B2}  Y(iJ) 

i=*l  u  u 

L  k3  k4K2 

+  l  1“  ?3  +  —  54)  l(lj) 


i  =  l 


K 


L  L 

x2i 


Y^k5 

+  l  6J*T(iL)  B5 }  Y(iL) 


K  L 


K 


Y2k5 


l  l  5?T(lj)  7  F(ij)  +  I  5**T(iL>  —  0(1) 

i=l  j»l  L  i=l  K 


(ID 


where,  as  it  can  be  seen  easily,  that 


1  1 

A  *  J  J  a  ^  aT  rr  d£dn 
0  0  ~ 


and 


*  "  {o  {o  ?T''>  dWn 


B3  =  /  a(0,n)aT(0,ri)dn  ,  B 2  =  1  a  r(0,n)aT  c(0,n)dn 
~  0  ~  ~  ~  0  ~  *  ~  9 


B3  =*  /  a(l,n)aT(l,n)dn  ,  84  =  /  a  f(l,n)aT  f(l,n)dn 

~  0  ~  ~  ~  0  ~  *  ~  y 

B5  -  J*  a(C  ,l)aT(5 ,0)dC 


F(ij)  -  J*  J*  a(C,n)6(ij)(C-i)dCdn  ,  G(1)  -  /*  a(£  ,l)Y(i)(Od£ 


(12) 
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Now  Eq.  (11)  can  be  assembled  in  a  global  matrix  equation 


6Y*t  K  Y  =  6Y*  F 


(13) 


By  virtue  of  the  fact  that  6Y*  is  not  subjected  to  any  constrained 


conditions,  one  has 


K  Y  =  F 


(14) 


which  can  be  solved  routinely.  Numerical  results  of  several  problems  in 
this  class  will  be  presented  in  a  later  section. 

III.  FORCE  VECTOR  DUE  TO  A  MOVING  CONCENTRATED  LOAD 

We  shall  describe  here  the  procedures  involved  to  arrive  at  the  force 
vector  contributed  by  a  moving  concentrated  load.  This  force  vector  has 
appeared  in  Eq.  (12)  as 


(15) 


The  shape  function  a(£,n)  is  a  vector  of  16  in  dimension.  In  the  present 
formulation  we  have  chosen  the  form: 


*  bi(5)bj(n)  ,  k  =  1,2,3, ...16 

if j  -  1,2, 3, 4 


(16) 


The  relations  between  k  and  i,j  are  given  in  Table  I.  These  are  the  conse¬ 
quences  of  the  choice  of  the  shape  function  such  that  Y(ij) ,  the  generalized 
coordinates  of  the  (ij)th  element,  represent  the  displacement,  slope,  veloc¬ 
ity,  and  angular  velocity  at  the  local  nodal  points.  Thus 


(17) 


P=1 

The  values  of  bjp  are  given  in  Table  II. 
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TABLE  I.  RELATIONSHIP  BETWEEN  (i,j)  AND  k  IN  EQUATION  (16) 


k 

1 

1 

- f 

(i,j)  1 

1 

k 

~1 

1 

1 

(i,j) 

1 

T- 

1 

1 

f 

(i,D  1 

9 

1 

1 

1 

(1,3) 

2 

1 

1 

1 

(2,1)  1 

| 

10 

1 

1 

1 

(2,3) 

3 

1 

1 

1 

(1,2)  1 
| 

11 

1 

1 

1 

(1,4) 

4 

1 

1 

| 

(2,2)  1 

| 

12 

1 

1 

1 

(2,4) 

5 

1 

1 

1 

(3,1)  1 

| 

13 

1 

1 

1 

(3,3) 

6 

1 

1 

| 

(4,1)  1 

| 

14 

1 

1 

1 

(4,3) 

7 

1 

1 

1 

(3,2)  1 

| 

15 

1 

1 

1 

(3,4) 

8 

1 

1 

1 

(4,2)  I 

1 

16 

1 

1 

1 

(4,4) 

TABLE  II.  VALUES  OF  t>ip  IN  EQUATION  (17) 


1 

1 

1 

P  ‘ 
i  1 

1  | 

2 

3 

1  4 

\ 

1 

1 

1 

1 

1  1 

1  I 

0 

-3 

|  1 

1 

1 

1 

1 

1 

| 

2  I 

0  1 

1 

-2 

|  1 

1 

1 

1 

1 

1 

1 

3  1 

0  1 

0 

3 

1  -2 

1 

1 

1 

1 

1 

1 

4  I 

0  1 

0 

-1 

|  1 

1 

1 

1 
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Now,  let  us  consider  ^(ij) (?"£)•  This  ’‘function"  represents  the 
effect  of  the  Dirac  delta  function  6 (x-x)  on  the  (ij)th  element.  If  the 
curve  of  travel  x  =  x(t)  does  not  go  through  the  element  (i,j),  ^ (± j ) ( 5—5 ) 

=  0.  If  It  passes  through  that  element,  one  has 

5(ij)(£  £)  =  <5(x-x)  =  K6(£-g)  (18a) 

with 

£  =  £(n)  (18b) 

The  function  £(n)  is  derived  from  x  =  x(t).  For  example,  if  the  force 
moves  with  a  constant  velocity,  one  has 

x  -  x(t)  =  vt  (19a) 

it  follows  from  Eqs.  (8)  that 

-  vK 

£  =  £(n)  =  -i+1  +  —  (n+j-1)  (19b) 

L 

With  Eqs.  (16),  (17),  (18),  and  (19),  one  writes  (15)  as 

,1  ,1  -  - 

F(ii)k  =  Kf  f  ak«fn)«(5-5)dWn  (20a) 

J  0  0 

.1  .1  -  -  p-1  q-1-  - 

F(ij)k  =  Kj  J  blpbjq  5  n  6(5-5)dMn 

00  (20b) 


Equation  (20)  can  then  be  evaluated  easily  once  the  exact  form  of  £  Is 
written.  For  example,  if  £  =  n ,  Eq.  (20)  reduces  to 
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F(ij)k 


4  4 

I  Ik  bjp  bjq 

p=l  q=l 


I 


1  ^p+q-2 

0 


dC 


4 

l 


q=l 


kbipbjq 

p+q-1 


(21) 


IV.  A  GUN  DYNAMICS  PROBLEM  AND  THE  MOVING  MASS  PROBLEM 

In  this  section,  the  solution  formulation  of  a  gun  tube  can  be  obtained 
as  a  special  case  to  the  gun  tube  motion  problem.  The  differential  equation 
of  this  problem  can  be  written  as. 

(Ely”)"  +  [P(x,t)y' ] '  +  pAy 
=  -  P(x,t)y”(x,t)H(x-x) 

-  mp[x2y”  +  2xy'  +  y]6(x-x) 

+  (rrip  g  cos  a)<5(x-x)  +  pA  g  cos  a  (22) 


The  notations  are  the  same  as  in  the  previous  section  if  they  have  already 
been  defined.  The  "gun  tube”  is  replacing  the  "beam"  whenever  appropriate. 

The  new  notations  are  defined  here: 

P(x,t)  =  ttR2 (x)p(t)  =  axial  force  in  the  tube  due  to  internal 
pressure  alone 

R(x)  *  inner  radius  of  tube 

p(t)  =  internal  pressure 

3j.  J.  Wu,  "A  Computer  Program  and  Approximate  Solution  Formulation  For  Gun 
Motions  Analysis,”  Technical  Report  ARLCB-TR-79019 ,  US  Army  Armament  Research 
and  Development  Command,  Benet  Weapons  Laboratory,  June  1979. 
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(23) 


,4  Jo  P^x 

P(x,t)  =  [ — P ( 0 ,t )  +  g( sin  a)  J  pAdx]  - 

°  pAdx 

o 

=  recoil  force  including  tube  inertia  in  axial  direction, 

H(x)  =  Heaviside  step  function 

x  =  x(t)  =  position  of  the  projectile 

mp  =  mass  of  projectile 

g  =  gravitational  acceleration 

a  =  angle  of  elevation 

With  similar  nondimensionalizat ion  as  before  and  assuming  that  the  cross- 
section  is  uniform,  ballistic  pressure  is  not  time  dependent.  Equation  can  be 
written  in  dimensionless  form 


yMH  +  [-P  +  g  sin  a] [ (l-x)yf ] f  +  Y2y 
=  -  P  y"H(x-x) 

-  Y2mp[x2y”  +  2xy?  +  y]6(x-x) 

+  mp  g(cos  a)6(x-x)  +  g(cos  a) 


(24) 


Where,  now,  everything  is  dimensionless  and 

2  C2  _  1  pA£4 
Y  T2  T2  El 


(25) 


It  is  also  clear  that  if  one  drops  the  second  term  on  the  left  hand  side  and 
the  first  and  the  last  terms  on  the  right  hand  side  in  the  above  equation,  the 
equation  becomes  that  for  a  moving  mass  problem. 
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A  variational  problem  associated  with  the  differential  equation  of 


Eq.  (24)  can  be  obtained  through  integration-by-parts. 


with 


12  3 

61  =  (6l)y  =  l  (6li)y  -  l  (SJj)  =  0  (26a) 

i-1  j-1 

A  A  -  A  A 

I 1  =  J  J  y"y*"dxdt  ;  Io  =  (P-g  sin  ct)  J  J  y'y*'dxdt 

0  0  z  0  0 

2  A  r1**  ~A  f1 

I3  =  -y/J  J  yy*dxdt  ;  I4  =  -PJ  J  y' y* ' H(x-x)dxdt 
0  0  0  0 

I5  =  -P/^  /^y'y*'6 (x-x)dxdt  ;  Ig  =  -mB2Y2/^  /^t2y'y*' 6(x-x)dxdt 
I7  =  -mg2Y2/^  f^ty'y* 1 6 ' (x-x)dxdt  ;  Ig  =  2m3 Y2/q  /  ty* y*6 (x-x)dxdt 

j  [ .  _  _ 

I9  =  -mY2/^  /^yy*6 (x-x)dxdt  ;  I^q  =  “niY2/^  /^yy*^ (x-x)dxdt 
I11  =  /  {k1y(0,t)y*(0,t)  +  k2y(0,t)y*(0,t)  + 


k3y( 1 >t)y*(l ,t)  +  k^y' (l,t)y*' (l,t) }dt 


I12  =  k7/Q  y(x,0)y*(x, l)dx 


(26b) 


and 


4  r1 

Ji  3  -g  cos  a  I  J  y*dxdt 
0  0 


A  A  -  - 

J2  *  -gm  cos  a  y*6(x-x)dxdt 


J3  =  ky/^  Y(x)y*(x,l)dx 


(26c) 
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The  variational  problem  also  produces  the  following  initial  and  boundary  value 
conditions  in  addition  to  the  differential  equation: 

y(x,0)  =  0 

•  -  1 

y(x,l)[l  +  m6(-  0-x)]  +  k.7[y(x,0)  -  Y(x)]  =  0  (27a) 

and 

y”(0,t)  -  k-2y'  (0,t)  =  0 
y"(l,t)  +  k4y' (l,t)  =  0 

y”'(0,t)  +  kiy(0,t)  +  (-P+g  cos  a)y'(0,t)  +  Py'(0,t)H(-  3t2) 

2 

+  meV (o,t)5(-  et2)  =  o 
2 

y” ' ( 1 >t)  -  k3(l,t)  +  Py'(l,t)H(-  et2-l)  +  mB2y'(l,t)6(-  0t2-l)  =  0 

2  2 

(27b) 

Other  than  the  fact  that  the  present  problem  is  much  more  complicated  than  the 
one  associated  with  a  moving  force,  the  basic  concept  of  solution  used  previ¬ 
ously  does  not  change  and  we  shall  omit  the  details  of  solution  formulation 
here. 

V.  NUMERICAL  DEMONSTRATIONS 

Some  numerical  results  obtained  will  now  be  presented.  Let  us  consider  a 
simply-supported  beam  subjected  to  a  unit  moving  force  with  a  constant 
velocity 

i 

v  =  — 

T 

As  T  varies  from  00  to  0,  the  velocity  varies  from  0  to  00 . 
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It  will  be  helpful  to  compare  v  with  some  reference  velocity  which  is  a 


characteristic  of  the  given  beam.  It  is  known  that  for  a  simply-supported 
beam,  the  first  mode  of  vibration  has  a  frequency  (see,  for  example, 


reference  4) 

w  1 

f1  =  —  =  —  ( — ) 

2ir  2ir  C 

TT  * 

=  —  (cycles  per  seconds) 

2C 

and  the  period, 

2C 

T],  =  “ 

TT 

where 

pA£4 

C2  - - 

El 

Consider  the  vibration  as  standing  waves.  They  travel  at  a  speed 

T\l 

vi  «  2*fi  =  — 

C 


Hence,  the  relative  velocity 

v  C  T1 

v  = 

vj  irT  2T 

IT 

We  shall  take  C  =  1.0  for  the  moving  force  problems.  Thus,  fj  ■  -  = 


1.5708  Hz.  T^  =  0.6366  sec.  and 

i 

ttT 

^K.  N.  Tong,  Theory  of  Mechanical 
257;  p.  308. 

Vibration,  John  Wiley,  New  York,  1960,  p. 

# 
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Using  a  grid  scheme  of  4x4,  Tables  III,  IV,  and  V  show  the  deflections 
as  the  concentrated  force  Q  =  1.0  moves  from  the  left  end  to  the  right  end  of 
the  beam.  Since  we  have  defined  that  T  is  the  time  required  for  the  load  to 
travel  from  one  end  to  another,  t  =  0.5T,  for  example,  denotes  the  point  when 
the  load  is  at  the  midspan  of  the  beam  if  v  is  constant.  At  t  =  1.0  T,  the 
force  has  reached  the  other  end  and  the  deflection  should  be  zero  in  the 
static  case. 

Solutions  by  Fouries  series  [1]  are  also  obtained  and  they  are  also  given 
in  these  tables  (numbers  in  parentheses)  for  close  comparisons. 

Table  III  shows  that  for  T  =  100  sec,  v  =  1/300  or  more  or  T  is  more  than 
300  times  the  natural  frequency  T^ ,  the  deflections  as  P  moves  across  the  beam 
is  nearly  the  static  deflection.  The  dynamic  effect  of  the  load  in  the  case  T 
=  100,  as  indicated  by  the  deflection  curve  at  t  -1.0  T  is  indiscernible.  For 
v  =  1/3  and  v  =  3.33,  the  dynamic  effect  is  very  much  pronounced  as  indicated 
by  Table  IV  and  V.  The  agreement  between  the  present  results  compared  reason¬ 
ably  well  with  the  series  solution  in  Tables  II  and  III.  It  is  extremely  well 
in  case  of  nearly  static  cases  as  shown  in  Table  III. 


15 


TABLE  III.  DEFLECTION  OF  A  SIMPLY  SUPPORTED  BEAM  UNDER  A  MOVING  LOAD 

(T  =*  100  sec.) 

y(x,t)/£  [x  10— 1 ] 


1  x/* 

It/T 

1 

1 

0. 

T 

1 

0.25 

1 

1 

0.50 

r 

i 

- r 

0.75 

r 

1 .00  I 

1 

1  o. 

1 

1 

0. 

T 

1 

0. 

1 

1 

0. 

i 

i 

~r 

o.  1 

“T 

o.  1 

1 

1 

1 

1 

(0.) 

1 

1 

(0.) 

1 

1 

(0.) 

i 

(0.)  1 

(0.)  1 

1 

I  0.25 

1 

1 

0. 

1 

1 

.1172 

1 

1 

.1432 

i 

i 

1 

.0911  I 

o.  1 

1 

| 

1 

1 

(0.) 

1 

1 

(.1167) 

1 

| 

(.1426) 

i 

i 

(.0907)  | 

j 

(0.)| 

1 

I  0.50 

1 

0. 

1 

1 

.1431 

1 

1 

.2082 

i 

i 

.1431  | 

o.  1 

1 

| 

1 

1 

(0.) 

1 

1 

(.1433) 

1 

| 

(.2085) 

i 

i 

(.1434)  | 

1 

(0.)l 

1 

|  0.75 

1 

1 

0. 

1 

1 

.0908 

1 

1 

.1427 

i 

i 

1 

.1168  I 

0.  1 

1 

1 

1 

1 

(0.) 

1 

1 

(.0916) 

1 

1 

(.1438) 

i 

i 

(.1176)  I 

1 

(0.)  1 

1 

1 

|  1.00 

1 

0. 

1 

1 

-.0047 

1 

1 

-.0066 

i 

i 

1 

-.0046  I 

o.  1 

1 

1 

1 

1 

(0.) 

1 

1 

(-.0002) 

1 

1 

(-.0003) 

i 

i 

(-.0002)  I 

1 

(0.)l 

1 

TABLE  IV. 

DEFLECTION 

OF 

A  SIMPLY  SUPPORTED  BEAM 

UNDER  A  MOVING  FORCE 

(T  =  1. 

.0 

sec . ) 

y(x,t)/£ 

[x 

10"1] 

1  x/£ 

i r 

T 

1 

i 

~r 

It/T 

i 

0. 

1 

0.25 

1 

0.50 

i 

0.75  I 

i  .oo  i 

1 

1  o. 

1 r 

i 

0. 

T 

1 

0. 

1 

1 

0. 

i 

i 

~r 

0.  1 

~T 

o.  1 

1 

I 

i 

l 

(0.) 

1 

| 

(0.) 

1 

1 

(0.) 

i 

i 

(0.)  1 

1 

(0.)  1 

1 

I  0.25 

i 

0. 

1 

1 

.09489 

1 

1 

.11349 

i 

i 

1 

.07108  | 

o.  1 

i 

i 

i 

1 

(0.) 

1 

1 

(.09795) 

1 

| 

(.11414) 

i 

i 

(.06942) I 

1 

(0.)l 

1 

|  0.50 

i 

0. 

1 

1 

.20542 

1 

1 

.30402 

i 

i 

1 

.21491  I 

o.  1 

i 

1 

i 

l 

(0.) 

1 

1 

(.20802) 

1 

1 

(.30257) 

i 

i 

(.21126) I 
j 

(0.)l 

1 

1 

|  0.75 

i 

0. 

1 

1 

.03869 

1 

.09522 

i 

i 

.09641  I 

o.  1 

i 

i 

i 

1 

(0.) 

1 

1 

(.05829) 

1 

1 

(.09397) 

i 

i 

(.08092) | 

j 

(0.)l 

1 

I  1.00 

i 

0. 

1 

-.10200 

1 

(.05191) 

i 

(.11574) | 

o.  1 

1 

1 

i 

i 

(0.) 

1 

1 

(.01994) 

1 

1 

(.03148) 

i 

i 

(.02405) I 

1 

(0.)l 

1 
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TABLE  V.  DEFLECTION  OF  A  SIMPLY  SUPPORTED  BEAM  UNDER  A  MOVING  FORCE 

(T  =  0.1  sec.) 

y(x,t )/l  [x  10"1] 


1  xM 

|  t/T 

1 

1 

0. 

1 

1 

0.25 

T 

1 

0.50 

T 

1 

r 

0.75  | 

1.001 

1  o. 

1 

1 

.0 

1 

1 

.0 

1 

1 

.0 

1 

1 

r 

•  0  I 

.0  I 

1 

(.0) 

1 

1 

(.0) 

1 

1 

(.0) 

1 

| 

(.0)  1 

1 

(.0)  1 

I  0.25 

1 

1 

.0 

1 

1 

.0619 

1 

1 

-.0148 

1 

1 

.0043  | 

•  0 

1 

1 

(.0) 

1 

1 

(.0645) 

1 

1 

(-.0149) 

1 

1 

(.0033)  | 

| 

(.0)1 

I  0.50 

1 

1 

.0 

1 

1 

.2002 

1 

1 

.1228 

1 

1 

-.0494  | 

•  0  1 

1 

1 

(.0) 

1 

1 

(.1952) 

1 

1 

(.1262) 

1 

| 

(-.0479)  | 

| 

(.0)1 

I  0.75 

1 

1 

.0 

1 

1 

.3007 

1 

1 

.3837 

1 

1 

.0770  | 

•  0  1 

1 

1 

(.0) 

1 

1 

(.2929) 

1 

1 

(.3849) 

1 

1 

(.0801)  j 

(.0)1 

1.00 

1 

1 

.0 

1 

1 

.4601 

1 

1 

.4912 

1 

1 

.5767  | 

.0  1 

1 

1 

(.0) 

1 

1 

(.4018) 

1 

1 

(.4880) 

1 

1 

(.5959)  | 

1 

(.0)1 
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FIGURE  1.  A  Typical  Finite  Element  Grid  Scheme  Showing 
the  (i,j)tk  Element  and  the  Global,  Local 
Coordinates. 
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